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A NEW DERIVATION OF THE EQUATIONS OF 
ALREADY-UNIFIED FIELD THEORY 


by 

Menahem Schiffer 
and 

Ronald Adler 

The algebraic Rainich conditions and the differential relations 
of Wheeler and Misner which form the basis of the already unified field 
theory are usually obtained from the Einstein-Maxvell system of equations 
by rather laborious tensor analysis. We here obtain these results by 
utilizing classical matrix theory and a special local coordinate system. 

0* Introduction 

The combined equations of classical vacuum electrodynamics and 
general relativity, 

*> V ■ * I * V 

(o.l) b) f^ v = 0 
llv 

c) *f^ v = 0 

llv 

are usually referred to as the Einstein-Maxvell equations. These were 
investigated by Rainich in 1925 [1] and Wheeler and Misner in 1957 [2,5], 
who found that the system (1) is equivalent to the following set of 
purely geometric relations: 



a) 


( 0 . 2 ) 


b) 

e) 

d) 


R 11 = 0 

V rV o ’ $ < R ,e rT@ ) <W 

R qo > 0 (system of real coordinates) 


v - v' = 0 ) 
X|t t|x 


vHIT 



Mv . 

1 _jt 

R** 


Subsequently, Wheeler, Misner and others used the system (2) as the 
basis of a geometrical-topological theory of gravitation and electro¬ 
magnetism, i.e. the already-unified field theory [2,5]. It is our 
purpose to give an alternative derivation of (2) from (1) which is 
based on classical matrix theory and i6 considerably simpler than the 
original derivation. 

1. Some properties of antisymmetric matrices 

T 

Let us consider a 4 X 4 antisymmetric matrix f = - f , with 
complex components. We will identify this matrix with the Minkowski 
electromagnetic field tensor in part 5. The characteristic polynomial 
0(X) = |XI - f| of f is easily shown to be an even function of X: 

(1.1) 0(X) « |XI - f| * |XI + f T | = |XI + f| = (-1) 4 |-XI - f| * 0(-X) 

Thus 0(X) = 0(-X), and it follows that 0 has the form 

(1-2) 0(X) = X 4 + a 2 X 2 + a Q 


2 



Furthermore, it is evident from (l.l) that if A, is an eigenvalue of f 
then -X is also. The eigenvalues, therefore, occur in two pairs: 

X^, -X v X 2 , -Xg. It i6 easy to express the coefficients and a Q 

in (1.2) in terms of these eigenvalues) since the eigenvalues are roots 
of 0(X), we can write 

( 1 . 3 ) 0(x) - (x-x 1 )(x+x 1 )(x-x 2 )(x+x 2 ) = \ k - (x 2 +x 2 ) x 2 +x 2 x 2 

Thus the coefficients a 2 and a^ are easily identified as 

(1-M » 2 a " < X 1 + ^ > a 0 x X 1 X 2 

The Cayley-Sylvester theorem tells us that f itself satisfies 
0(f) » 0, so we have from (1.2) and (1.4) 

(1.5) t k + a 2 f 2 + a Q = t k - (X 2 + X 2 ) f 2 + X 2 X 2 = 0 

The preceding results allow ub to construct an interesting 
symmetric matrix from f. Define 

(1-6) T = f 2 + -j- I = f 2 - \ (X 2 + X 2 ) I 

This will be identified in part 3 with the energy-momentum tensor of 
the electromagnetic field. From the definition (1.6) and equation (1.5) 
it follows immediately that T 2 ^Ls. a multiple of the identity matrix : 
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(1-7) 


T 2 


1 (X 2 
k (X 1 



I 


2 2 p 

We will always assume that X^ ^ X^, so that T is not zero. 

A second interesting property of T follows from a consideration 
of its Jordan canonical form for similarity. By Jordan's theorem any 
complex matrix is similar to a direct sum matrix as follows; 


(1.8) T = Q _1 ^Q = i 

where the have scalars 

first superdiagonal: 



0 



on the diagonal and "l"s on the 



However, for the present case the are severely limited by the 

requirement (1.7). Indeed, if one squares the equation (1.8) and 

compares with (1.7) it is evident that the must all be 1 x 1 

12 2 

matrices and that the must have the values _+ — (X^ * Xg). 
Thus 


k 



( 1 . 10 ) 


T 


-1 



Q 


and the are clearly the eigenvalues of T. 

The signs of the can be determined from the definition 

(1.6) and the fact that they are the eigenvalues of T: 


(1.11) | T i I - T| = |[T i - ± (X 2 ♦ X 2 ))I - f^| = 0 

12 2 p 

Thus t - — (X^ + Xg) is an eigenvalue of r and must consequently 
2 2 

be X^ or Xg, each occurring twice. This gives 


( 1 - 12 ) 


*1 = Tg = - T = - = g (X^ ” Xg) - T 


and upon substitution into (1.10) 


(1.15) T - Q -1 ^Q - Q -1 [ 1 ° ) Q ; 

V 0 -Tl / 


I = 


l o 


0 1 


Thus it is evident that T has a null trace and can be written as 


(1.14) 


- f 2 - r (Tr f 2 ): 


5 



In part 4 we will need the fact that the Q which appears in 

"1 T 

(1.13) may he chosen to be orthogonal: Q = Q . The proof of this is 
straightforward. Since T is symmetric, (1.13) gives 

(1.15) T T = Q T ^(q‘ 1 ) T = T = Q" 1 ^ Q 


from which 


( 1 . 16 ) f/(QQ T ) = (Q Q T ) //' 

T 

If QQ is now written in terms of 2X2 submatrices as 


(1-17) 



7 

0 


then substitution into (1.16) reveals that y = 6 = 0. Thus 


(1-18) 



and the submatrices a and 0 are clearly symmetric and have nonzero 
determinants. 

Next note that the choice of Q is somewhat arbitrary since 
it may be replaced by 


(1.19) 


Q » 



Q. 
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where R and S are arbitrary 2X2 non6ingular matrices. This is 
evident from the fact that 


( 1 . 20 ) 


T = Q = Q* 1 


Furthermore, it is evident that 


( 1 . 21 ) 



Thus in order to complete the proof we need only 6how that R and S 
may be chosen so that R0( R^" = SO S T = I, for then by (1.20) and (1.21) 
will be the desired orthogonal matrix. This is an easy task} let 


( 1 - 22 ) 




T 

and substitute into R3 R * I. This gives three equations in the four 
unknowns u, v, w, and z. 


2 2 
au + 2buv + cv =1 

(1.25) aw 2 + 2bwz + cz 2 = 1 

auw + buz + bwv + czv = 0 


If a is nonzero, a solution is 


( 1 . 24 ) 




w = 


-b 


^afaT 
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The case of nonzero c is completely analogous. If both a and e 
however, are zero we can use the following solution: 


(1.24 1 ) 


z = 1 


2b 


v = - i 


w = 


JL_ 

2b 


This exhausts all cases, so we have shown that R does exist. In 

completely similar fashion there exists a nonsingular S such that 

T 

Sp 3 =1, so the proof of the statement is complete: the Q in 

(1.13) can be chosen to be orthogonal, 

The canonical form (1.15) for T can be utilized to invests 
gate the structure of its generating matrix f. Write f in the form 


(1.25) 



Q 


where Q i6 the same orthogonal matrix as in (1.13) and K, L, M, 
and N are to be determined Since f is antisymmetric, both K and 
N are also antisymmetric. From the definition of T (1.6) and the 
canonical form (1.13), one finds easily that 


( 1 , 26 ) 


f 2 ,q t 



Q 


Using (1.25) for f and (1.26) for f 2 and the obvious identity 
^f - ff 2 = 0, we obtain 


(1 27) 


(x 2 - X 2 )L = (X 2 - x 2 ) M « 0 


8 



2 2 

Since X^ and X 2 are assumed to be unequal, L and M are zero. 
Similarly, the identity ff ■ f 2 yields 

(1.28) K 2 = X 2 I ; N 2 = X 2 I 


Since K and N are antisymmetric this tells us that 


(1.29) 


K = iX x J ; 


N = iX 2 J ; 



1 

0 


) 


The signs in (1.29) are arbitrary and have been chosen to be positive. 
This gives the following canonical form for f 

/iX^J 0 

(1.5°) t - « T ( o ^ 

where Q is an orthogonal matrix. 

2. A simplification of preceding results 

Having obtained the canonical form of f in (1.30) we can 
put all the preceding results in very simple form. Define 



( 2 . 1 ) 





Q 


where Q is the same orthogonal matrix as in (1.30). Then one finds 
by an elementary calculation that 


9 



( 2 . 2 ) 


2 2 3 3 

p + q - i ; p = p ; q'= q 

4 k 

qp = pq ■ 0 j p + q ■ I etc. 

The canonical form (1.30) nov is expressible as 

(2.3) t = X L p + l 2 q 

From the definition of T and (2.2), we have 

(2 A) T - \ (X* - X 2 )(p 2 - q 2 ) 

It immediately follows that T _is_ traceless and that T? is a 
multiple of the Identity 

(2-5) T* - £ U 2 - X 2 )V + q 4 ) * £ (Tr T 2 ) I 

These are the principal results of part 1, but now made quite trans¬ 
parent . 

The degree of uniqueness in the relation of f to T is 
quite interesting and easily ascertained. If we construct a T matrix 
fran a different ? defined by 

(2.6) f = X^p + X 2 q j X^ = ^X 2 - X 2 cosh a 

X 2 = A 2 - X 2 sinh a 


10 



where a ie an arbitrary parameter, then 

(2.7) ? - £ (X 2 - X 2 )(p 2 - q 2 ) - | (X 2 - X 2 )(p 2 - q 2 ) - T 

that is, f and f generate the same T independent of the choice 
of a. Note in particular that the choice 

^1 ^2 

(2.8) cosh a * ■ —- ; einh a = — - 

A 2 - x| /x 2 - x| 

clearly yields the original f, i.e. ? = f. It is, therefore, clear 
that an entire one-parameter family of f matrices generates the same 
T matrix. 

3- The algebraic Rainlch conditions 

We now wish to apply the results of parts 1 and 2 to the 
task of deriving the Rainich conditions (0.2a, b, c). For convenience 
we will work in a locally geodesic system, so that at some fixed world- 
point P the Christoffel symbols vanish and the ordinary and co¬ 
variant derivatives of first order coincide. Such a locally geodesic 
system is only determined up to a linear transformation with constant 
coefficients. Thus, in order to use matrix theory, we may use a 

geodesic system in which the metric tensor g is the Kronecker S 

° uv 

This allows us to ignore the distinction between contravariant and 
covariant indices and makes tensor algebra and matrix algebra the 
same. Two further features of this special system should be noted; 


11 




firstly, it is only unique up to an orthogonal transformation, which 
leaves 8 invariant. Secondly, in such a system both the coordinates 
x^ and the Minkowski tensor f will in general be complex. 

With the above choice of a coordinate system, we can write the 
energy-momentum tensor of the electromagnetic field (0.1a) in matrix 
notation as 

(3.1) T « f 2 - £ (Tr f 2 ) I 

Comparison with (1.16) showB that T depends on f in precisely the 
same way as the T matrix which we investigated in part 1 depended on 
the antisymmetric matrix f. Thus all the results concerning T in 
part 1 are immediately applicable to the electromagnetic energy- 
momentum tensor. In particular, we can assert that T is traceless 
and T 2 is a multiple of I. By equation (0.1a) we can say the same 
about R . That is, in tensor notation, 

(5.S) W 

Furthermore, since (3-2) is written in covariant form it is true at all 

world-points and in all coordinate systems. 

It is possible to strengthen equation (3-2b) by demonstrating 

that, under reasonable physical assumptions (as explained below), the 
t6 

scalar R is positive and real. This is easily shown as follows. 

From (2.5) we know that R x p R T ^ may be expressed in terms of the 

2 2 2 

eigenvalues of f^ as (X^ - X g ) . This is a covariant statement. 
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Indeed, the eigenvalue equation of may he written covariantly as 

(3.5) - w a - s p 

from which it is clear that the eigenvalue X is a scalar. Equation 
(3*3) gives rise in the usual way to a covariant secular equation for X: 


( 5.10 


' f ap ‘ X| 


= o 


What we wish to 6how now is that if we make the physically reasonable 

demand that f^ be real in a system of real coordinates then 
2 2 2 

(Xi - X 2 ) is a positive real scalar . 

A scalar can be calculated in any coordinate system so we will 
momentarily utilize a real tangent Lorentz system with 


(3-5) 




-1 


We can write f^, which is now assumed real, in the general form 


(3-6) 


ap 
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and for convenience we introduce the real three-tuple vectors 
C = (a, b, c) and D = (f, -e, d). The secular equation (3*4) with 

from (3-5) and f^ from (3-6) is found by direct expansion to be 

(3.7) X 4 - (C 2 - eP) X 2 - (C-D) 2 - 0 

where we have used three-dimensional vector notation. This equation 
has the two solutions 

(3.8) X 2 = C 2 - D 2 ± /C 2 - D 2 ) 2 * 4(C-P) 2 
from which we find immediately that 

(3-9) R tP R tP - U 2 - x|) 2 * (c 2 - rP) 2 + U(C-D) 2 

2 2 

By our standing assumption that X^ X^ this is obviously a positive 
real number. Since it is also a scalar, we have shown that it is real 
and positive in general. 

The above fact is interesting as a statement concerning (3*2b), 
but it is also important in the derivation of condition (2.b) as we 
shall see presently. 

It is a well-known fact that in special relativity the component 
T^ of the energy-momentum tensor is proportional to the energy density 
of the electromagnetic field and must, therefore, be a positive number. 
In order to preserve the Bame interpretation of T qo in general rela¬ 
tivity we must show that 
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(3.10) 


R 


oo 


> 0 


is a consistent covariant demand, at least in all real coordinate 
systems. 

B 

To demonstrate the consistency of this relation consider R Q 
as a linear transformation on an arbitrary vector v R 


(3.11) 


V B 
Ot p 


From (3.2b) it is evident that 


(3.12) 


¥<",15 


r tP ) 


Since R^ R T ^ ha6 been shown to be positive, we see that w Q and v Q 
are both timelike, both spacelike, or both null. In particular, we 
can say that R a P carries the light cone into itself. 

Next we consider R in a new coordinate system 


(3.13) 


dx 4 dx V R 

as 0 &° 


Let be v^ and R 

as 0 as 0 ^ 

be written 


be w as above, so that R can 
p 7 oo 


(3-14) 
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Thus it is clear that the demand that S be positive definite is 

oo 

equivalent to the covariant requirement that carry a vector v Q 

into a vector w° such that v°w > 0 • This can be interpreted 
physically by 6aying that carries the forward light cone into 

itself and the backward light cone into itself. Thu6, by the statement 
(3.10), one gives a sign sense to in a covariant way for any real 

coordinate system. 

The conditions (3.2) and (3-10) were first obtained by Rainich 
in 1925> using somewhat different matrix methods than we have used. 

They are usually referred to as the algebraic Rainich conditions and, 
indeed, are algebraic in the sense that they follow from a consideration 
of f and T at a single world-point P. 

uv nv 

This completes the derivation of (0.2a, b, c) from (0.1), so 
only the differential relation (0.2d) remains to be considered. This 
differential relation will be our most important and interesting result. 
We have derived the algebraic conditions (0.2a, b, c) mainly for the 
sake of completeness and in order to lay a groundwork for the final 
relation (0.2d). 

4. The differential relations of Wheeler and Mlsner 

In order to obtain (1.2d) we need to take into account the fact 
that f obeys the Maxwell equations (0.1b,c). At the world-point 
P in our special system we can ignore the distinction between contra- 
variant and covariant indices and between ordinary and covariant 
derivatives of first order. This allows us to rewrite (0.1b,c) as 


16 



(4.1) a) f = 0 j b) *f - 0 

nv|v M.V I V 

where we have retained the Einstein summation convention without regard 
for index position. 

Our aim now is to study how equations (4.1) reflect themselves in 
properties of the tensor R . Observe first that, given an energy mo- 
memtum tensor according to (2.6), one can only obtain up to 

a parameter a which may depend on the point considered. Algebraically, 
the ct field at different world-points could be completely incoherent, 
but we will find that the Maxwell equations (4.1) provide a differential 
system for a which makes it a determined point function. The inte- 
grability condition on the system will be the Wheeler-Misner relation 
which we seek. 

Let r * A^ - Xg and write f in the general form (2.6). 

(4.2) f = r cosh a p + r Binh a g 

4V UV UV 

By inspection of (2.1) the dual of f is found to be 

(4.3) #f * r cosh a g + r sinh a p 

(JLV uv uv 

Recall that our special coordinate system is still arbitrary to 
an orthogonal transformation. Thi6 allows us to assume without loss of 
generality that the matrix Q which appears in the definition of p 
and q is precisely the identity matrix I. Then p and q at the 
world-point P take the particularly simple forms 
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(4.4) 



It follows then from (2.4) that the energy-momentum tensor becomes 


(4.5) 



and that the Maxwell equations at P are 


(4.6) 


a) f = p (r cosh a) + q (sinh a) 
4v|v |v Iv 


+ r cosh Op + r sinh a q * 0 
4v|v nv|v 


b) *f - q (r cosh a) + p (r sinh a) 
nv|v ^ |v ^ 


+ (r cosh a) q + (r sinh a) p 

4 v|v nv|v 


0 


To simplify (4.6) define at the world-point P vectors "Tf 
and K and scalars A and B as follows: 

n 


(4.7) TT = P ; K=q j A = r cosh a j B = r sinh a 
4 4 v|v nv|v 


Then Maxwell's equations (4.6) reduce to 


(4.8) 


a) p A + q B 

nv | y 


Iv 


b) 


q, „ A + p B 

V , v V J, 


+ A TT + B K 

I 'n \x 

+ A K + B TT 
4 V 


= 0 

= 0 
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Now multiply a) by p, b) by q, add the two results, and use (2.2) 
to obtain easily 

(M> Atp^ TT V + K v + (106 r) ] + Bt\ v TT/ P„ v \ + o ! - 0 

II ^ 

Similarly, multiply a) by q, b) by p, and add to get 

(4.10) A[^ v TT V + p^ v K v + o ] + B[^ y TT V + K v + (log r)^) * 0 

By the definition of A and B and the structure of these equations, 
it is clear that (4.9) and (4.10) are consistent only if the coeffi¬ 
cients of A and B are zero. Thus we obtain 


a) 


(4.11) 


a * - (q TT + p k ) 
I v %iv v v 


b) (log r) = - (P uv TF V + V K v ) 


As we indicated above the a field is now seen to be subjected to a 
simple differential system. 

The next step in our development is to calculate the vectors 
TT^ and and 4o use the result to put (4.11) into a more inform¬ 
ative and interesting form. In order to do this we first will calculate 
the derivatives of the matrices p and q. At the world-point P 
these matrices are precisely the p and q in (4.4), but as we move 
to a nearby world-point P(e;n), displaced from P by a distance € 


in the (i-th coordinate direction, these matrices will became p 


M 
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and q^ which differ from p and q by a small rotation in space- 
time corresponding to an orthogonal matrix . Let us represent 


this matrix by a series 

(*• 12 ) 


U) * 1 * 8 c (n) * 8 V) + 


^ j oucii uu^iicD uiau ^(ll) 

symmetric, as is well known, and that its inverse is 


(*•13) 


) * ^"i) = 1 ■ € c (^) + ■■■ 


It is now easy to calculate the derivatives of p and q at 
P in the u~th direction. We have from (2.1), (4.12), and (*.13) 

l 4 - 14 ) » (ll) - <4) »«<„)**♦ *l» ' %) » 1 + °(' 2 > 

Now subtract p from both sides, divide by e, and pass to the limit 
e -» 0 to obtain 

(>.•15) f |u - »<=(,)- V) » 

and similarly for q . If we write in terms of 2x2 sub¬ 

matrices as 


(*•16) 
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then the antisymmetry of implies that ajT^ = - a^, 

T T 

d^ • - d^, and b^ = - c^. Now substitute (4.16) and (4.4) 
into (4.15) to obtain an explicit form for p 


f 0 ^(n) 


p. = i 


'* V“ c (n) J 0 


Similarly, we find for q 


V' 1 


0 -V) J 


JC M 0 


We will later need the derivative T at P, so we will 

In 

obtain it now while it is most convenient. Recall that T may be 


written as 


m 1 2, 2 2. 

T = o r (P - q ) 


Using (4.4), (4.5), (4.17), and (4.18) it is then easily seen that 


T = (log r) 2 T + r 

In In 


0 


Observe that in the above expressions for p , q , and T 

In In In 
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only the 2X2 matrices 
write these explicitly as 


occur. Let us 
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A straightforward calculation then yields 



We have achieved an explicit form for K and TT in terms of 

n 'n 

the components of . Let us substitute now the results (4.25) into 

(4.11a) and use (4.4) to get an explicit form for a 

In 


(4.24) 




Observe an important fact at this point: the vector a 


is composed of the same components k^, m^, and n^ which 

occur in T and is, therefore, expressible in terms of T alone. 

In 

Indeed, our only remaining problem is to express the correspondence of 


a and T in covariant form. To do this consider the covariant 

In In 

vector 
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Ot.25) 


l 




T T 0 * 

px 


which we will now calculate at P in our special coordinate system. 

By (4.5) the denominator is simply r\ The numerator takes the simple 
form 


(4.26) 




MV y 

T T ' 


e T 

xv&r 


T „ 
04|v 07 


Using the explicit form (4.20) for T we have 

Mv 


(4.27) * 2(log r) 


XVP7 04|V * 7 


X.V07 04 47 


0 b, V 

2 M S 

r € XV07 . -C^) 0 T 47 


yve 


The Rainich condition (3.2b) and the Einstein equations (0.1a) tell us 

that T^ T is a multiple of the metric tensor and is, therefore, 

symmetric in 0 and 7 ; since e.. _ is totally antisymmetric the 

A,VP7 

first term of (4.28) vanishes identically. A short calculation using 
(4.21) and (4.5) then reveals that 


(4.28) 
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and a slightly tedious but elementary calculation gives 


(4.29) 



Comparing this with (4.24) ve see that 


(4.30) 



__ Mv 


T 


7 

Ji 


which is in tensor form and therefore valid in general. Using the 
Einstein equations (0.1a), we can also write this in geometric form as 


(4.31) 






rNI 


V 


U 


Pk 


The differential relations of Wheeler and Misner follow 
immediately from (U. 31 ), for in order that (4.31) be an integrable 
relation we must have 

(4.32) v - V * 0 

x|t t|x 

Thus we have finally obtained (0.2d), the last basic relation of the 
already-unified field theory. 
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